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We demonstrate spectroscopy and thermometry of individual motional modes in a mesoscopic 2D 
ion array using entanglement-induced decoherence as a method of transduction. Our system is a 
^400 /xm-diameter planar crystal of several hundred ^Be^ ions exhibiting complex drumhead modes 
in the confining potential of a Penning trap. Exploiting precise control over the ^Be^ valence electron 
spins, we apply a homogeneous spin-dependent optical dipole force to excite arbitrary transverse 
modes with an effective wavelength approaching the interparticle spacing (^-^20 /xm). Center-of-mass 
displacements below 1 nm are detected via entanglement of spin and motional degrees of freedom. 
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Studies of quantum physics at the interface of micro- 
scopic and mesoscopic regimes have recently focused on 
the observation of quantum coherent phenomena in op- 
tomechanical systems [1-3 . The realization of quan- 
tum coherence in mechanical oscillations involving many 
particles behaving approximately as a continuum pro- 
vides exciting insights into the quantum-classical tran- 
sition. Previous work has shown that crystals of cold, 
trapped ions behave as atomic-scale nanomechanical os- 
cillators [4]-[6], with the benefits of in-situ tunable mo- 
tional modes and exploitable single-particle quantum de- 
grees of freedom (e.g. valence electron spin). Our system 
of hundreds of crystallized ions in a Penning trap provides 
a bottom-up approach to studying mesoscopic quantum 
coherence. In this context, the relevant particle numbers 
are sufficiently small to permit excellent quantum control 
without sacrificing continuum mechanical features. Be- 
yond these capabilites, trapped ions have long provided 
a laboratory platform for studying diverse physical phe- 
nomena including: strongly-coupled one-component plas- 
mas (OCPs) [IllH]; quantum computation [9l[T0] and sim- 
ulation pTHT5] : dynamical decoupling [16]; and atomic 
clocks and precision measurement [17 . 

In this Letter, we present an experimental and theoret- 
ical study of motional drumhead modes in a 2D crystal 
of ^Be+ ions confined within a Penning trap. We ex- 
cite inhomogeneous modes of arbitrary wavelength (see 
Fig. [ija)) through application of a homogeneous, spin- 
state-dependent optical dipole force (ODF) to a large- 
scale spin superposition. Distinct drumhead modes are 
entangled with the ^Be+ valence electron spins by tun- 
ing a beat frequency {fiR) between two ODF lasers near 
a mode resonance. This spin- motion entanglement is de- 
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tected as a /ii^-dependent decoherence of ion spins whose 
magnitude conveys the specific mode temperature. 




FIG. 1: (color online) (a) Calculated structure of selected 
transverse eigenmodes (bm) for a 2D crystal of 331 ^Be^ 
ions. Mode frequencies, cjm, decrease as effective wave- 
length gets shorter. The arbitrary color scale indicates rel- 
ative ion displacement amplitude. One example of an ion 
spin state with similar symmetry is given below each of the 
four highest-frequency eigenmodes. The symbol x (•) denotes 
spin-projection into (out of) the plane. Interaction between 
these spin and mode configurations mediated by the spin- 
dependent optical dipole force (ODF) leads to excitation of 
the corresponding eigenmode. (b) Illustration of a single 
plane of ^Be^ within the Penning trap. Two 313-nm beams 
intersect at the ion cloud to form a traveling wave of beat 
frequency fiR and effective wavevector along the direc- 
tion of the trap magnetic field. The electric field intensity 
is uniform in the plane, but the spin-dependent induced AC 
Stark shift permits excitation of transverse modes of arbitrary 
wavelength. 
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Previous global mode studies on 2D planar ion ar- 
rays were restricted to modes with wavelengths on the 
order of the cloud size p!8H22] . By contrast, the short- 
wavelength modes studied here are of particular interest 
due to their increased sensitivity to strong-correlation 
corrections [23l [24] compared to those with long wave- 
length, which are well-described by fluid theory. Ther- 
mometry of large Coulomb crystals has thus far been 
limited to determination of global temperature through 
Doppler profile measurements [25 , which give a mini- 
mum sensitivity of ~0.5 mK in ^Be+. Our temperature 
measurement is mode-specific and may be employed be- 
low the Doppler cooling limit, providing an alternative 
to Raman sideband thermometry [26 . 

The Penning trap used for this work is detailed in 
a previous publication [27|. Application of static volt- 
ages to a stack of cylindrical electrodes provides har- 
monic confinement along z (the trap symmetry axis) 
with a ^Be+ center-of-mass (COM) oscillation frequency 
of uj\l2i\ = 795 kHz that is independent of the num- 
ber of trapped ions. The trap resides within the room- 
temperature bore of a superconducting magnet, and ra- 
dial confinement is achieved via the Lorentz force gener- 
ated by rotation of the ion cloud through the static, ho- 
mogeneous magnetic {B) field of ~4.46 T oriented along 
z. Application of a time-dependent quadrupole 'rotating 
wall' potential permits phase-stable control of the rota- 
tion frequency (cJr), and thus the confining radial force of 
the trap [28l [29] . In the limit of a weak rotating wall po- 
tential, the harmonic trap potential in a frame rotating 
at Ur is [8] 

9$trap(r,z) = ^Mco! {z^ + I3r^) , (1) 



where M (q) is the mass (charge) of a single ^Be+, 
= 27r X 7.6 MHz is the cyclotron frequency, and z 
(r) is axial (radial) distance from the trap center. We set 
the rotation frequency, cj^, such that the radial confine- 
ment is weak relative to transverse confinement (/3 <C 1), 
resulting in a single ion plane. 

The mj = ±1/2 projections of the Be+ ^5'i/2 ground 
state are split by ~124 GHz and serve as 1 1) and ||) 
'qubit' states, respectively. Global spin rotations are per- 
formed by injecting 124-GHz radiation through a waveg- 
uide attached to the side of the trap. The ^Be+ ions 
are Doppler laser cooled with laser beams directed both 
parallel and perpendicular to z. Both beams are tuned 
to the ^6'i/2(mj = +l/2)-^P3/2(mj = +3/2) transition 
at ~313 nm to cool ion motion below 1 mK. This same 
transition is used for ion detection and projective spin- 
state measurement. Discrimination of |t) (bright) from 
\i) (dark) is performed with a fidelity > 99% [27 . 

The axial and radial confining potentials are tuned 
to yield a planar ion configuration. Due to mutual 
Coulomb repulsion and the low ion temperature, the ions' 
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FIG. 2: (color online) (a) Pulse sequence used for excitation 
and detection of transverse motional modes. Global spin rota- 
tions are performed with microwaves at ^-^124 GHz, while the 
state-dependent optical dipole force is applied in each arm of 
the spin echo for a duration r. We implement 7r-pulse times 
(tvr) as short as 65 /iS. (b) Measured (points with statistical 
error bars) and fit (solid blue line) probability of detecting 
It) (^t) of the spin echo sequence. Frequency- 

dependent decoherence is due to entanglement of spins with 
the axial COM mode {uji/27t = 795 kHz) as a function of 
ODF detuning Si = /j.r — uji in a cloud of 190 ± 8 ions. Each 
point is an average of 90 experimental runs. The fit provides 
a mode temperature of 2.3 di 0.5 mK, whose error includes a 
5% uncertainty in ODF beam angle, fe. For comparison, the 
lower (upper) dashed line is calculated assuming 0.4 mK (4.0 
mK). (c) Illustrated phase-space trajectories of state W)n at 
different detunings, in a frame rotating at cji. Axis la- 
bels represent COM momentum {pz oc Im[aji]) and position 
{z oc Re[aji]). 



minimum-energy configuration is a 2D crystal with trian- 
gular order [30 . Ion spacing is ~20 /im, and individual 
ions can be resolved using stroboscopic imaging at uOj.- 
The planar array of TV ions exhibits 3 A/" motional modes, 
N of which are drumhead oscillations transverse to the 
crystal plane (see Fig. [ija)). As with ID ion strings, 
the frequencies of these transverse modes decrease with 
decreasing effective wavelength due to screening of con- 
fining electric fields by nearby ions. The transverse eigen- 
vectors ihm^ ^ ^ [I5 ^]) corresponding eigenfrequen- 
cies iujrn) are obtained by first numerically calculating 
the zero-temperature 2D ion configuration in the pres- 
ence of the Penning trap potentials. Applying a Taylor 
expansion of the combined trap and Coulomb potential 
about each ion equilibrium position, we diagonalize the 
N X N stiffness matrix whose eigenvalues and unit eigen- 
vectors are Um and bm , respectively [311 [32] . The rela- 
tive displacement amplitude of an ion j is given by the 
jth element of 6^, denoted as where l^jmf — 

To excite transverse modes in our 2D Coulomb crystal, 
we employ a spin-dependent ODF generated by interfer- 
ing two off-resonant laser beams at the ion cloud posi- 
tion. This is depicted schematically in Fig. [ijb). The 
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two ODF beams are produced from a single beam us- 
ing a 50/50 beamsplitter and subsequently pass through 
separate acousto-optic modulators that allow fast (~1 
/is) switching and impart a relative detuning juji. The 
beams intersect at an angle of 6>r = 4.8° ±0.2° at the ion 
cloud position, and their relative alignment is adjusted 
to orient the effective wavevector (A^) of the resulting 
standing {/j^r = 0) or traveling {/j^r ^ 0) wave to within 
~ 0.05° of z. The common wavelength (313.133 nm) and 
unique linear polarizations of the beams are chosen such 
that the AC Stark shift from the interfering beams on 
state It) is equal in magnitude and opposite in sign to 
that on \D [33 . The result of the interference between 
these two beams is a spin-dependent force on each ion, j 
{F^j = —F^j = Fj). The Hamiltonian for this interac- 
tion is ^ODF = — Y^^=i FjZj{t) COS {fiRt) a j ^ where Zj{t) 
is the time-dependent position operator and (7| is the z- 
component Pauli operator for ion j [14]. The elliptical 
beam waists (100 /im x 1000 /im, with the major axis ori- 
ented parallel to the ion plane [33]) are sufficiently large 
to generate an approximately uniform ODF with varia- 
tion below 10% across the ~400 /im-diameter planar ion 
crystal. Typical ODFs for this work are Fj ^ 10~^^ N 
along z. 

Figure [2|a) illustrates the experimental control se- 
quence for microwaves (black line) and ODF lasers 
(shaded regions) used to coherently excite transverse 
modes of motion. Ions are ffist prepared in the 'bright' 
state It) AT = Y[f=i Itj) via optical pumping [27]. The se- 
quence of microwave pulses in Fig. [2|a) comprises a spin 
echo (SE) [34 that, in the absence of the ODF beams, 
rotates the ions to the 'dark' state \i)N with >99% fi- 
delity. The SE cancels low-frequency precession about 
z due to ODF laser intensity and magnetic field fluctu- 
ations as well as microwave phase noise [^ The 
spin-dependent ODF is applied in each arm of the SE for 
a duration r. 



The initial microwave pulse rotates each spin by 7r/2 
to produce the state YljLi (Itj) ~ which is a su- 
perposition of all possible (2^) spin permutations. Im- 
portantly, it is the creation of this state that permits 
subsequent excitation of arbitrary transverse modes with 
our homogeneous, spin-dependent ODF. By tuning jur 
near a mode of frequency cj^, the spin-dependent ODF 
excites those components of the spin superposition with 
approximately the same symmetry as the eigenvector bm- 
A subset of these eigenvectors and associated spin states 
are illustrated in Fig. [TJa). Depending on experimental 
parameters, the spin states may be entangled with differ- 
ent motional states at the end of the control sequence of 
Fig.[2|a). Upon measurement of the spin state (perform- 
ing a trace over the motion), entanglement is manifested 
as spin decoherence that varies with fiR. We observe this 
as a decrease in the length of the spins' Bloch vector and 
a concomitant increase in the probability (P^) of measur- 
ing state It) averaged over all ions. 

Figure [2|b) gives experimental and theoretical results 
for a sweep of fiR near the COM frequency, cji, with 
r = 500 /is and Si = {fiR — ui). On resonance (^i = 0), 
the pulse sequence leads to excitation (de-excitation) of 
the COM mode in the ffist (second) arm. When the 
product |(5ir/27r| = Hs a non-zero integer, each spin state 
traverses / full loops in phase space over r (see Fig.|2|c)). 
At intermediate detunings, the spin and motion remain 
entangled at the end of the pulse sequence, producing 
the lineshape of Fig. |2|b). These motional excitations 
are described by the spin-dependent displacement opera- 
tor U{t) = n,,mexp [{ajmal - a]^am)^]\ [3lE3[36], 
where ot^rnij) is the coherently-driven complex displace- 
ment amplitude for ion j of mode m, and aj^(a^) is the 
creation (annihilation) operator for mode m. Accounting 
for both arms of the pulse sequence, we obtain [33] 



^ — [c^m(l - cos(/)) + iiJ,Rsin(f)- e*"^"^^ {um [cos (fiRr) - cos {ij,rt + (p)] - ijiR [sin {iirt) - sin {iirt + </))]}] , 

(3) 



where h is Planck's constant, Fj is the ODF magnitude 
on ion j, and (/) = (r + tj^){jj.R — ujm) accounts for phase 
evolution of the ODF drive relative to that of the mode. 

Although the coherently driven, spin-dependent dis- 
placements (o^jm) are independent of the initial motional 
state (assuming Lamb-Dicke confinement [14 ), the spin- 
motion entanglement signal in Fig. pjb) sensitively de- 
pends on this initial state. This can loe qualitatively un- 
derstood in terms of the spatial structure of a harmonic 
oscillator Fock state, |nm), of mode m. A state |n^) ex- 
hibits n wavefunction nodes and therefore, as n increases, 
a fixed spin-dependent displacement results in less wave- 
function overlap between different spin components due 



to the increasing spatial frequency of |n^) wavefunctions. 
This leads to larger decoherence and greater displace- 
ment sensitivity as the average mode occupation, n^, 
is increased for a given mode. We ffi the experimental 
measurements in Fig. |2|b) using theory that attributes 
a thermal state of motion to each mode m characterized 
by mode occupation ~ kBTm{fi^m)~^ and tempera- 
ture T^. Neglecting spin-spin correlation contributions, 

we find the probability P^^"^ of detecting ion j in state |t) 

at the end of the pulse sequence to be [33] 



exp -2^\ajm\'^{2nm + 1) 



(4) 
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FIG. 3: (color online) (a) Measured (lower) and calculated (offset) probabilities for measuring |t) after the spin echo sequence 
as a function of ODF beat frequency for a sweep of /j^r over the first five transverse modes with 250 ±15 ions. The modes at (jJ2 
and cjs are split due to distortion of the ion cloud boundary by the rotating wall potential. Panels (b) and (c) give results of 
wider sweeps with cjr/Svr = 43.2 kHz and 44.7 kHz, respectively, in a crystal of 345 ± 25 ions. Frequency-dependent deviation 
from ^0.1 is due to spin-motional entanglement, while the background is due to spontaneous emission from the ODF beams. 
The histogram (red bars) shown below each experimental curve depicts the density of calculated eigenmodes at the given ujr. 
Histogram bins are 10 kHz wide and plotted with an arbitrary vertical scale. As described in Fig. T]^a), the highest-frequency 
feature is that of the COM mode and the ^50 lowest- frequency eigenmodes include nearest-neighbor ions oscillating out of 
phase. Features at cjr and precise harmonics thereof (shaded in light green) are due to spin-motion entanglement with in-plane 
degrees of freedom excited by the small (^ 10~^ Fj) component of ODF perpendicular to z. 



Here F accounts for decoherence due to spontaneous 
emission induced by the ODF lasers over the duration 2r, 
and is responsible for the background level of P-^ ~ 0.1 ob- 
served in all experimental data presented here [37 . The 
total detection probability is obtained by averaging 

all P^^\ 

For interaction with the COM mode {bji = ;^,Vj G 

[1,A^]), aji is obtained from Eq. ^ through measure- 
ment of the ODF laser intensities ^4j and trapped-ion 
number, while F is determined from decoherence ob- 
served with /j^R detuned far from any modes. As such, 
the only parameter of Eq. Q not measured directly is ni, 
which is varied to fit experimental data as in Fig. [2|b) , 
where we obtain ni = 60 ± 13 (Ti = 2.3 ± 0.5 mK). 

We note that a detectable phase-space displacement is 
obtained with a very small amplitude of |Qfjrn|- For ex- 
ample, in Fig.|2|b), the 20% decrease in the Bloch vector 
at ^ir/27r ^ ±1.4 corresponds to a spin-state-dependent 
excitation of the COM mode with a mean excursion of 
~0.6 nm in each arm of the pulse sequence. This shift is 
less than 0.2% of the wavefunction spread of a single ion 
in the planar array. Our sensitivity to displacements im- 
proves with increasing mode temperature provided that 
the ODF is adjusted to avoid full decoherence {P^ = 0.5) 
at the detuning of interest. 

Figure [sj^a) shows the result of a sweep of /hr over five 
transverse modes and corresponding theory. The theo- 
retical spectrum (offset for clarity) is generated assum- 
ing Ti = 10 mK and Tm>i = 0.4 mK, with Ti obtained 
from a fit. The large COM temperature of Fig. [sja) is 
produced by quickly switching off the i-oriented Doppler 
cooling beam on a time scale of ~27rcj^^. In this case, 
sudden loss of radiation pressure from the cooling light 
induces a COM oscillation amplitude of ~50 nm that we 
detect as an elevated ni. A more adiabatic reduction of 



the cooling beam intensity yields ni ~ 26 (Ti ~ 1 mK). 
For modes other than the COM, we must additionally 
calculate the bjm values for the trap potentials and ion 
number in a given experiment. For these modes, we find 
temperatures consistent with the Doppler cooling limit 
of 0.43 mK. 

To measure the full spectrum of transverse modes, we 
repeat the sequence of Fig. [2|a) for 30 kHz < /j.r /27r < 
800 kHz with r = 1 ms. With the exception of the 
COM mode, the frequencies of the remaining N — 1 
modes depend sensitively on our choice of crystal rota- 
tion frequency, Ur [20 . Figures [3]^b)-(c) show the result 
of these experimental runs for uJr/27T = 43.2 kHz and 
44.7 kHz, respectively. For this ion number of 345 ± 25, 
the single-plane configuration is stable over the range 
42.2 kHz < ujr/27r < 45.2 kHz. Histograms of calcu- 
lated mode density versus (ir/ 27: are plotted below each 
experimental curve with an arbitrary vertical scale and 
bin width of 10 kHz. The distribution of eigenfrequen- 
cies narrows as uOj. is decreased; weaker radial confine- 
ment (see Eq. Q) leads to lower ion densities and re- 
duced screening of trap potentials, thereby moving the 
frequency of the shortest-wavelength mode toward that 
of the COM. This behavior is clearly visible in Figs.jsj^b)- 
(c). Additionally, we find quantitative agreement be- 
tween the measured spectrum and that generated from 
numerical calculation of the transverse eigenmodes under 
the given experimental conditions, documenting coupling 
to both short- and long-wavelength modes. The sharp 
features of Figs. |3jb)-(c) shaded in light green reflect ex- 
citation of in-plane resonances at harmonics of uJr due 
to a very small component of the ODF (~10~^Fj) along 
the ion plane. These spectral features may be reduced 
through more careful alignment of to i, but their 
strong response suggests an elevated motional tempera- 
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ture perpendicular to z. 

In summary, we have used entanglement of spin and 
motional degrees of freedom to map the full transverse 
mode spectrum of a mesoscopic 2D ion array. This 
technique provides a tool for sensitively and accurately 
measuring the temperature and displacement amplitude 
of individual drumhead modes, facilitating identification 
of mode-specific heating mechanisms and the resulting 
non-equilibrium energy distributions. Coherent, spin- 
dependent excitation of transverse modes is the basis for 
engineering quantum spin-spin interactions with trapped 
ions pTI4T4[ [3T |l32|l38| l39]. making mode characterization 
a critical element of such experiments. Future work will 
include investigation of low-frequency in-plane modes at 
frequencies smaller than uOj.- A predicted subset of these 



modes includes in-plane shearing motion whose restoring 
force is due exclusively to strong correlations. 
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Supplementary Material: 
Spectroscopy and Thermometry of Drumhead Modes in a Mesoscopic Trapped-Ion Crystal using 

Entanglement 



Optical Dipole Force Details 

Figure [4] shows a simple sketch of the optical dipole force (ODF) laser beam set-up. As discussed below, the 
frequency as well as the beam polarizations were chosen to null the AC Stark shift from an individual beam and 
to produce a state- dependent force which is equal in magnitude but opposite in sign for the |t) and |^) qubit states 
{F^ = —F^ ). The off-resonant laser beam frequency was detuned from the cycling transition {^3/2^ = 3/2)) 

by Ar —63.8 GHz. This gives detunings of 15.6 GHz and 26.1 GHz respectively from the |t) |^3/27^J = 1/2) 
and 1^) 1^3/27 = ~l/2) transitions. Laser beam waists were Wz — 100 /im in the vertical (z-direction) and 
Wx — 1 mm in the horizontal direction. Here we define the waist as the distance from the center of the beam over 
which the electric field intensity decreases by 1/e^ (i.e. I{z) ~ e~^^/^^^ ). With the small 2.4^ incident angle each 
beam makes with respect to the plane of the crystal, this provided greater than 90% uniform electric field intensity 
across ion crystal arrays with N < 250. 

We used linearly polarized laser beams. Let 

Fu {r, t) euFu cos {ku • r- uut^ 
EL{r,t) clFl cos (^kL ' r- ULt J 

denote the electric fields of the upper and lower ODF beams. If (f)p is the angle of the laser beam electric-field polar- 
ization with respect to vertical polarization (e • x = 0) , then the AC Stark shift of the qubit states when illuminated 
by a single beam can be written 

At, acss = cos^ i^p) + sin^ i^p) 
A;, acss = cos^ ((/>p) + sin^ 

where A^{A^) is the Stark shift of the |t)(|i)) state for a 7r-polarized beam (e parallel to the z-axis) and B^{B^) is 
the Stark shift of the |t)(|i)) state for a a-polarized beam (e perpendicular to the i-axis). (Here we neglect the small 
a polarization (oc sin(2.4^)) that exists when (j)p = 0.) The Stark shift of the qubit transition is 



{A^ - A^) cos^ (^p) + {B^ - B^) sin^ (0^) 



(5) 



a) 



kjj =k\ 0-cos — ,-sin^ 






kr =k\ 0,- COS — , sin 



b) 



FIG. 4: Sketch of ODF laser beam setup, a) The ODF laser beams lie in the y-z plane at angles zLOr/2 with respect to the 
?/-axis. b) View looking in the —y direction. The beams are linearly polarized but with different polarization angles relative to 
vertical polarization. 
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FIG. 5: Sketch of the ID optical lattice wave fronts (red lines) generated by the ODF laser beams. These wave fronts need to 
be aligned with with the ion planar array (represented by the blue dots). Here Xr = 27r/|A^| ^ 3.7 /xm and Oerr denotes the 
angle of misalignment, do ^ 400 fim is the typical array diameter for N ~ 200 ions. With the wave front alignment technique 
discussed in the text we obtain ^err < 0.05°. 



If — and — B^ have opposite signs, there is an angle which makes Aq 



0. For a laser detuning of 



-63.8 GHz, A, 



at 



±65^ 



With 

^acss — for each ODF laser beam, we exploit the freedom to choose their polarization in order to obtain 
a state-dependent force. Specifically we choose Eu to have a polarization given by (j)p^u = 65^ and El to have a 
polarization given by (j)p^i = —65^. In this case the interference term in the expression for the electric field intensity 

(^Eu + El) produces a polarization gradient which results in spatially dependent AC Stark shifts 

(A^ cos^ {(j)p) — B^ sin^ {4^p)) 2 sin {5k • z + /URt) 
{A^ cos^ {(j)p) - B^ sin^ {(l)p)) 2 sin {5k ■ z ^ iirI) 



for the qubit levels. Here 5k = 

beams, /j^r = uu — ujl is the ODF beat note, and (j)p 
produces a state-dependent force F^^i{z^t) = F^^^ cos {5k 



2 A: sin (^) is the wave vector difference between the two ODF laser 

|0p,n| = The spatially dependent AC Stark shift 

— fiRt) where 



F^ = 2 5k {A^ cos^ {(l)p) - B^ sin^ ((/)p)) 
F^ = 2 5k {A^ cos^ {c^p) - B^ sin^ {c^p)) ' 



F^. We operate at A^^ = —63.8 GHz where for Aacss = we also obtain F^ 
^p,u 1 4^p,i 1 and A^^ we use straight forward atomic physics along with well known values for the energy 



In general F-^ ^ 
For a given 

levels and matrix elements of ^Be+ to calculate F as a function of the electric field intensity Ir 
at the center of the laser beams. For Or = 4.8^ and Ir = 1 W/cm^ , F = 1.5 x 10~^^ N. 



\Eu\ 



Wave Front Alignment 

The ODF laser beams produce a ID optical lattice characterized by the effective wave vector 5k and beat note 
liiR. In the previous section we assumed that 5k \\ z, or equivalent ly that the wave fronts of the lattice were aligned 
perpendicular to the i-axis (magnetic field axis). If the wave fronts are not normal to the z-axis as sketched in Fig. [sj 
then the time dependence of the optical dipole force seen by an ion in the rotating frame depends on the (x, y) position 
of the ion. This complicates the interaction generated by the optical dipole force and is avoided by careful alignment. 

We used top-view images (images of the ion resonance fluorescence scattered along the magnetic field) from a 
single plane to measure a misalignment of the ODF wave fronts. For this measurement we set /ir = (stationary ID 
lattice) and detune the frequency of the ODF laser beams approximately 0.5 GHz below the |t) |^^3/2 = +3/2) 
Doppler cooling transition. This small detuning generates sufficiently large AC Stark shifts on the cooling transition 
to measurably change the ion scatter rate from the Doppler cooling laser. With the Doppler cooling laser on and the 
ODF beams turned off we observe a spatially uniform, time-averaged image of a rotating planar crystal. With the 
ODF beams on, ions located in regions of high electric field intensity at the anti-nodes of the optical lattice are Stark 
shifted out of resonance with the Doppler cooling laser. This is what produced the dark bands in the top-view image 



FIG. 6: Top- view image of the spatially inhomogeneous fluorescence from a single ion plane produced by the AC Stark from a 
static (/jR = 0) optical dipole force lattice with misaligned wave fronts. Dark bands are regions of high standing wave electric 
field intensity (parallel to the dashed yellow line) . The bright horizontal feature bisecting the center of the image is fluorescence 
from the weak Doppler laser cooling beam directed perpendicular to the magnetic field. The image was obtained by subtracting 
a background image with the ODF beams off. 



shown in Fig. [6] From images like this we determine how to move the ODF beams to align the wave fronts normal to 
z. Improved alignment is indicated by a longer wavelength fringe pattern. With this technique we have aligned the 
ODF wave fronts with the planar array to better than O^rr ^ 0.05^. 

Images like that shown in Fig. |6] were typically obtained with 1 s integration. This means the imprint of the ID 
lattice on the planar arrays was stable during the integration time and indicates a phase stability of our ID lattice 
of better than Is. We note that direct fluorescence imaging of the ID lattice, for example by tuning the ODF laser 
resonant with the Doppler cooling transition, is not viable. Even at low powers, resonantly scattered photons across 
the large horizontal waist of the ODF beams apply a large torque, causing the rotation frequency and radius of the 
array to rapidly change, typically driving the ions into very large radial orbits. 



Spin-Motion Entanglement Produced by the Spin-Dependent Optical Dipole Force 

With the wave vector 6k of the ID optical lattice aligned parallel to the optical dipole force generated by the 
lattice is independent of the ion position and can be written 

F^{t) = -F^{t)=F cos i^Rt) (6) 

where jur is the frequency difference between the ODF laser beams. More generally we allow for the possibility that 
the ODF laser intensity could be diferent for each ion, resulting in a different spin-dependent force Fj for each ion j, 

F,^{t) = -F,^{t) = F, cos {^iRt) . (7) 

In the experimental set-up, the variation in Fj is less than 20%. The ODF interaction with the ion spins can be 
written as 

N 

HoDF = -^Fj COS {iiRt) Zja] . (8) 
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Here zj is the axial position operator for the j ion, which can be written in terms of the axial normal modes 
^6rn,^m^ of the planar array, 



N 



— ^ 



m—1 



2Mujr, 



^ ale'''-') 



(9) 



The eigenvectors are normalized so that \^jm\^ = Xlj l^jml^ = 1- Both the eigenvectors hm and eigenfrequencies 
ujrn are calculated by solving for the ion equilibrium positions and diagonalizing the stiffness matrix obtained by 
Taylor expansion of the potential about the ion equilibrium positions [1]. 

The Hamiltonian Hqdf of Eq. ([8| is time dependent. The evolution operator for Hqdf is obtained from a second 
order expansion of the Magnus formula [U |3] 



UoDF {t) = exp 



HoDF{t')dt' 



2f9 



t rt 
dt2 



rt2 

Jo 



ODF {h) , Hqdf {ti)] dti 



(10) 



Higher order terms do not contribute as the commutator [HoDF{'t2), HoDF{ti)] commutes with HoDF{t'). Fohowing 
the discussion of Ref. [3_, UoDF{t) can be written 



UoDF{t) = exp 
= exp 



Ei (Em (ajm(i)aL - o:*^ii)am) <t|) 
UsM{t)-Uss{t) 



■ exp 



(11) 



The first term UsM{t) describes spin-dependent displacements ajm{t) of the normal modes m where, for the cos(/ii^t) 
time dependence of the interaction in Eq. ([8|, 



[uJrn - e'"^""^ (^m cos(/ii^t) - i/J.R sm{/j.Rt))] . 



(12) 



Here zom = \/fi/{2Mujm) - The second term Uss{t) describes an effective spin-spin interaction where the pairwise 
coupling Jj^k{t) is given by 



FjFk ^ bjmbkmZom j sin(/ii^ - UJm)t ^ OJm sin(/ii^ + UJrn)t _ UJm sin(2/ij^t) 
2^2 ,,2_ _„,2 



^^R 



I^R - ^r) 



(13) 



For now we assume Uss{t) can be neglected. We will discuss the validity of this assumption at the end of this section. 

The interaction UsM{t) = exp (X]m {^jm{t)dl^ — Q^j^(^)^m) generates spin- motion entanglement that is 
the subject of this study. The commutator 





because is real. Therefore we can write UsM{t) as a product of individual spin displacements 



(14) 



which is Eq. (4) of the Letter. By neglecting the spin-spin entanglement {Uss{t)) we can independently calculate the 
evolution of each spin j. 

We now calculate the spin motion entanglement generated by Usnif) during the free precession period of a Ramsey 
sequence shown in Fig.[7|^a). The calculation for the spin-echo sequence of Fig.jTj^b) used in the experiments is identical 
except for a more complicated expression for the aj^(t)'s (see next section)?^ Each spin j is prepared in state W) 
at the start of the sequence. If an ODF is not applied during the free precession period, the spin is rotated to the 
dark \\) state by the final 7r/2 pulse of the sequence. With the application of a spin-dependent ODF, in general the 
spin is entangled with the motion at the end of the Ramsey sequence. We detect this spin-motion entanglement by 
measuring the probability of finding spin j in the |t) state. Let 



^5m(^) = exp ^ {ajm{t)al - a]^{t)am) 
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a) 



71/2 




ODF(hr) 



71/2 



''^II) 


n/2\ 


y ^ X 




71/2 1 








ODF(^r) J 


ODF(^r) 










X 


X 







FIG. 7: Pulse sequences described in this supplemental material, a) Ramsey pulse sequence consisting of two 7r/2 rotations 
with an intermediate arm of duration r during which the ODF is applied, b) The spin echo sequence repeated from Fig. 2(a) 
of the Letter which consists of two arms of duration r. 



denote the evolution of spin j by the spin-dependent ODF. By re- writing 



^5m(^) = cosh {ajmCii, - aj^cimij + sinh (a^ma^ - ^jm^m)j cr] 



we calculate 



sinh (ajm^L - ^jm^m)^ ^ siuh (ajm^L " ^jrn^m)^ 



(15) 



th 



where Pi^sM denotes the probability of measuring the |t) state for spin j produced by the UsM{t) interaction, and 



( denotes an expectation value averaged over a thermal (Maxwell-Boltzmann) distribution of modes. We evaluate 
Eq. (15) by writing the sinh functions in exponential form. It is then necessary to evaluate expressions of the form 

e^e^) where A and B are operators which are linear in the raising and lowering operators aj^ and am- In this 

/ th 

case we can make use of the result 



/ th 



to obtain 



1 - exp ( -2^\ajm(t)\^ {2nr, 



1) 



(16) 



Here kBTm/ {fujJm) is the mean occupation number of a Maxwell-Boltzmann distribution characterized by 

temperature T^. We measure the probability of detecting |t) averaged over all the ions - 

The simple result of Eq. (16) was obtained under the assumption that we could neglect Uss{t) in Eq. (11). In 

general Uss{t) will contribute to the measured P^^\ This can be straight forwardly estimated when /ir is tuned 
close to the COM mode uji. In this case the resulting pair- wise interaction coefficients are identical for all ion pairs 
JjAt) - J{t) with 



J{t) 



F'2 



2n? N 



^01 



cji sin(/ii^ — co'i)^ cji sin(/ii^ + co'i)^ ujisin{2jj.iit) 



■ UJit 



For small detunings \fiR — \ <C cJi, J{t) is approximately bounded by \J{t) \ ^ J -t where 



J 



^01 



The fully connected, uniform Ising interaction exp 



iJ 



obtained by coupling through the COM mode 



is identical to the single-axis twisting interaction analyzed by Kitagawa and Ueda [5|. We use the expressions given 
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in Ref. [5] to calculate P^^ls^ probability of measuring spin j in the |t) state at the end of the Ramsey sequence 
due to the Uss{t) interaction, 



[NS{Jt)^ 



(17) 



This expression is valid for short times t where P^^lg is small. 

We obtain strong spin- motion entanglement for small detunings \fiR — uJi \ <C cji. The magnitude of the coherently 
driven amplitude aj,m=i(^) in the expression for UsM{t) (Eq. (14)) and P^^Im (^^- (1^)) maximized for a detuning 

— uji \ Tv/t where 



a, I t 



ll^R -^l| 



Fzoi 



-2uji 



The above expression neglects terms of order {juji — ojm) /^m- Inserting 
exponent is small gives 



into Eq. (16) and assuming the 



D(i) 



2KilL. (2ni + l)) 



(18) 



We compare P^^l^ (Eq. ^) with P^^lj^ (Eq. (18)) 



pU) 
p(j) 



N • 8 {jty 



^01 



2KilL.(2ni+l) 4fi2 2ni + l 



(19) 



For the work reported here F ^ 10 N, zqi = Y^^^TI^McJiy ~ 30 nm, and ni ~ 10 (Doppler cooling limit). For 
a typical interaction time t < 10~ s we calculate 'P^ s s ^ '^'\ s ivi 

< 0.1. Therefore for small detunings satisfying 
I Mi? — ^1 1 ^ (27r)/t <C uji we expect the spin-motion entanglement signature generated by UsM{t) to dominate contri- 
butions due to Uss{t)' We note that the spin- motion entanglement signature (^P^^sm^ decreases with temperature. 
For ground state cooling it may not be possible to neglect Us sit)- 

We do not estimate P^^ls for fiR tuned close to modes ujm other than the COM mode uji. Therefore we do not 
know if it is a good approximation to neglect Uss{t) when resonantly coupling to non-COM modes. However, we 
experimentally observe that neglecting Uss{t) gives a good description of our experimental measurements for /j^r 
tuned close to the tilt {uj2 and us) and the next lower frequency modes {u^ and u^). 



Spin Echo Sequence with Decoherence 

To calculate ajm{t) for the full spin echo sequence used in the experiment (see Fig. mb)), we must account for the 
accumulated phase difference between the ODF drive and oscillating ion cloud over the first arm and intermediate 
microwave 7r-pulse of combined duration (r + 1^^). This requires derivation of ajm{t) for an ODF interaction with an 
arbitrary phase offset, (/), given by the more general 

AT 

HodfW = -^Pj cos {^Rt + (t)) zfG], (20) 



where = (r + t-K^l^R — ^m) = (t + t-K^^m- Following the previous derivation of ajm{t) for = (Eq. (12)), we 
obtain 

(^jm{t, (/>) = ^^^-^^^Q^ [^^ cos((/)) - ijiiR sin(^) - e*"^^^ {cjm cos(/ii^t + (/>)- i/HR sin(/ii^t + (/))}] . (21) 



We will now define a new that may be substituted for ajm in Eq. (14) to calculate P^ for the full spin echo 
sequence exhibiting arm durations of r: 



<^fm = Q^jm(r, (/> = 0) - ajm{r, ^), (22) 
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where the above expression is given exphcitly in Eq. (6) of the Letter. 

To justify implementation of Eq. (22), it is useful to calculate for a single spin undergoing both the Ramsey 
and spin echo sequences. To simplify notation, we define the displacement operator D{ajm) = exp(aj^aj^ — aj^dm) 
which is applied separately to (g) |?/^^) and (g) IV^m), where is an arbitrary motional state of mode m. 



Assuming the state is initialized to g) IV^m), we calculate the result of the Ramsey sequence, P. 



(j) Ramsey 



to be 



p(j) Ramsey _ 

- 4 



[d^ (aj^ir, 4>)) - {-ajm{T,4>))} {h.c.} 



(23) 



where {h.c.} denotes the Hermitian conjugate of the first bracketed expression. Here the arbitrary phase (j) has no 
physical significance since its value is common to all displacements, and we have once again made the assumption 
that Fj^ = —Fj^. However, the spin echo result given by Pi"^^^^ is 



{d^ {-ajm{r, (/))) Z)t {ajm{r, 0)) - (a,^(r, (/>)) (-a,^(r, 0))} {h.c.} \^m) (24) 
\{^^^\ {Z)t (af^) - Z)t (-af^) } {h.c.} \^^). (25) 



We obtain Eq. (26) from Eq. (25) using the multiplicative properties of D and neglecting overall phase factors that 
leave Pi"^^^^ unchanged. Note that Eq. (26) is identical to Eq. (24) after an appropriate redefinition of ajm. 



Finally, the derivation of Eq. (16) neglected the effects of spontaneous emission from the ODF laser beams. De- 



coherence of the Bloch vector due to spontaneous emission from off-resonant light is well studied in our system [6]. 
The qubit levels are closed under spontaneous light scattering; that is, spontaneous light scattering does not optically 
pump an ion to a different ground state level outside of the two qubit levels. In the presence of off-resonant laser 
light, the decrease in the Bloch vector due to spontaneous scattering during the arms of a spin-echo sequence is 



']',spon 2 '■ 



1 



■ exp (-r • 2r)] 



Here F = {TRam + Tg^) /2 has contributions from both Raman scattering and elastic Rayleigh scattering that can be 
calculated from the laser beam parameters. With the spin echo sequence, we account for spontaneous emission by 
modifying Eq. ( 16 ) as follows 



)(i)SE 



1 - e-^^^ exp -2 ^ I' (2n^ + 1) 



^ t,SM 2 

where r is the length of time of a single arm of the spin-echo sequence. 
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